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ABSTRACT 

We present detailed integral field unit (IFU) observations of the central few kilopar- 
secs of the ionised nebulae surrounding four active central cluster galaxies (CCGs) 
in cooling flow clusters (Abell 0496, 0780, 1644 and 2052). Our sample consists of 
CCGs with Ha filaments, and have existing data from the X-ray regime available. 
Here, we present the detailed optical emission-line (and simultaneous absorption line) 
data over a broad wavelength range to probe the dominant ionisation processes, exci- 
tation sources, morphology and kinematics of the hot gas (as well as the morphology 
and kinematics of the stars). This, combined with the other multiwavelength data, will 
form a complete view of the different phases (hot and cold gas and stars) and how they 
interact in the processes of star formation and feedback detected in central galaxies in 
cooling flow clusters, as well as the influence of the host cluster. We derive the optical 
dust extinction maps of the four nebulae. We also derive a range of different kinematic 
properties, given the small sample size. For Abell 0496 and 0780, we find that the stars 
and gas are kinematically decoupled, and in the case of Abell 1644 we find that these 
components are aligned. For Abell 2052, we find that the gaseous components show 
rotation even though no rotation is apparent in the stellar components. To the degree 
that our spatial resolution reveals, it appears that all the optical forbidden and hydro- 
gen recombination lines originate in the same gas for all the galaxies. Based on optical 
diagnostic ratios ([OHI]A5007/H/3 against [NH]A6584/Ha, [Sn]AA6717,6731/Ha, and 
[OI]A6300/Ha), all galaxies show extended LINER emission, but that at least one 
has significant Seyfert emission areas, and at least one other has significant HH like 
emission line ratios for many pixels. We also show that the hardness of the ionising 
continuum do not decrease with galactocentric distance within our field-of-view (with 
the exception of one galaxy that show a core separation) as the emission line ratios 
do not vary significantly with radius. This also indicates that the derived nebular 
properties are spatially homogeneous. We fit AGN and pAGB stars photoionisation 
models as well as shock excitation models to our derived diagnostic ratios. These fits, 
combined with information from multiwavelength studies reveal that AGN photoioni- 
sation is the most likely ionisation mechanism in at least two cases even though shocks 
and pAGB stars can not be conclusively eliminated. 

Key words: galaxies: formation - galaxies: elliptical and lenticular, cD - galax- 
ies:clusters:individual:Abell 0496 - galaxies:clusters:individual:Abell 0780 - galax- 
ies: clusters: individual: Abell 1644 - galaxies: clusters: individual: Abell 2052 



1 INTRODUCTION 

Decades ago, elliptical galaxies were thought to contain very 
little, if any, gas. Studies of galaxy formation, therefore, of- 
ten focussed on the stellar properties, however we now know 
that a large fraction of the baryonic mass in massive galax- 
ies is believe to be in diffuse form. Thus a complete view 
of galaxy formation and evolution necessarily incorporates 
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both the stars and hot gas and an understanding of the 
processes by which these phases interact (McCarthy et al. 
2010). 

Cooling-flow clusters are common in the local Universe 
and massive central cluster galaxies (CCGs) are often found 
at the centres of these systems (Peres et al. 1998). If the 
central cluster density is high enough, intracluster gas can 
condense and form stars at the bottom of the potential well. 
Since the radiative cooling times for intracluster gas are 
short enough that gas can cool and settle to the cluster cen- 
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tre (Edge, Stewart & Fabian 1992), it has been suggested 
that the big envelopes of CCGs may arise from the grad- 
ual deposition of this cool gas. More recently, high spectral 
resolution XMM-Newton observations showed that the X- 
ray gas in cluster centres does not cool significantly below a 
threshold temperature of kT ~ 1 — 2 keV (Jordan et al. 2004, 
and references therein) . This initially contradicted the model 
that these young stars are formed in cooling flows. However, 
it is possible that star formation is ongoing in cool-core clus- 
ters at a much reduced rate (Bildfell et al. 2008). 

Previous studies have reported several examples of on- 
going star formation in CCGs, in particular those hosted by 
cooling-flow clusters (Cardiel, Gorgas & Aragon- Salamanca 
1998; Crawford et al. 1999; McNamara et al. 2006; Edwards 
et al. 2007; O'Dea et al. 2008; Bildfefl et al. 2008; Pipino et 
al. 2009; Loubser et al. 2009). However, the origin of the gas 
fuelling this star formation is not yet known. Possible expla- 
nations include processes involving cooling flows or cold gas 
deposited during a merging event (Bildfell et al. 2008). These 
processes will leave different imprints in the dynamical prop- 
erties, the detailed chemical abundances, and the star forma- 
tion histories of these galaxies, which can be studied using 
high-quality spectroscopy (Loubser et al. 2008; 2009; 2012). 

Observations that support this idea are blue- and 
ultraviolet-colour (UV-colour) excesses observed (indicative 
of star formation) in the central galaxy of Abell 1795 by Mc- 
Namara et al. (1996) and molecular gas detected in 10 out of 
32 central cluster galaxies by Salome Sz Combes (2003). The 
observations by Cardiel et al. (1998) were consistent with 
an evolutionary sequence in which star formation bursts, 
triggered by radio sources, take place several times during 
the lifetime of the cooling flow in the centre of the cluster. 
However, McNamara & O'Connell (1992) found only colour 
anomalies with small amplitudes, implying star formation 
rates that account for at most a few percent of the mate- 
rial that is cooling and accreting onto the central galaxy. 
Cooling-flow models for CCG formation also imply the for- 
mation of larger numbers of new stars, for which there is 
no good observational evidence (Athanassoula, Garijo & 
Garcia-Gomez 2001). The central cluster galaxies often host 
radio-loud AGN which may account for the necessary heat- 
ing to counteract radiative cooling (Von der Linden et al. 
2007). 

In summary, cooling flow models predict more cooled 
gas than is observed (Bohringer et al. 2001). Thus, it is pos- 
sible that the mass deposited into the molecular clouds is 
heated by one of several processes - hot young stellar pop- 
ulations, radio- loud AGN, X-rays or heat conduction from 
the intracluster medium itself, shocks and turbulent mixing 
layers, and cosmic rays. Therefore, only a small fraction of 
the cooled gas is detected (Crawford et al. 2005; Ferland et 
al. 2009). Thus, CCGs lie at the interface where it is cru- 
cial to understand the role of feedback and accretion in star 
formation. Within these cooling-flow CCGs, cool molecular 
clouds, warm ionized hydrogen, and the cooling intracluster 
medium are related. A complete view of the star formation 
process incorporates the stars with the gas and an under- 
standing of the processes by which these phases interact, 
and therefore, requires information from several wavelength 
regimes. 

In conclusion, although CCGs are probably not com- 
pletely formed in cooling flows, the flows play an important 



role in regulating the rate at which gas cools at the centres 
of groups and clusters. In the ACDM cosmology it is now 
understood that local massive haloes assemble late through 
the merging of smaller systems. In this picture, cooling flows 
seem to be the main fuel for galaxy mass growth at high red- 
shift. This source is removed only at low redshifts in group 
or cluster environments, due to AGN feedback (De Lucia 
&; Blaizot 2007; Voit et al. 2008). Indeed, if AGN feedback 
is not properly assumed in hydrodynamical simulations, an 
apparently bluer BCG is formed as a result of an accelerated 
late stellar birthrate even after the epoch of quiescent star 
formation (Romeo et al. 2008). 

We proposed and obtained IFU observations of the cen- 
tral few kiloparsecs of the ionised nebulae in active CCGs in 
cooling-flow clusters. These observations will map the mor- 
phology, kinematics and ionisation state of the nebulae to 
gain an understanding of their formation, heating and re- 
lationship to the cluster centre. We selected galaxies from 
the McDonald et al. (2010) study, which consisted of an Ha 
survey of 23 cooling flow clusters. Amongst their conclu- 
sions, McDonald, Veilleux & Rupke (2012) find a stong cor- 
relation between the Ha luminosity contained in filaments 
and the X-ray cooling flow rate of the cluster, suggesting 
that the warm, ionised gas is linked to the cooling flow. We 
chose objects with confirmed extended Ha emission, and 
with near-IR (2MASS), ultraviolet (GALEX), X-ray data 
(Chandra), and in some cases VLA 1.4 GHz fiuxes, already 
available (McDonald et al. 2010). Detailed properties of the 
host clusters, which are reported to influence the activity in 
the central galaxy, such as central cooling times and the off- 
set between the cluster X-ray peak and the central galaxy, 
have been derived. 

Parage et al. (2010) presented IFU observations of one 
nearby BCG showing LINER-like emission, and Brough et 
al. (2011) presented IFU data on four CCGs at z ~ 0.1 to cal- 
culate the dynamical masses of CCGs and measure their stel- 
lar angular momentum (although none of their four CCGs 
contained emission lines). The SAURON sample includes 
only one CCG (M87). Hatch, Crawford & Fabian (2007) 
and Edwards et al. (2009) analysed data comparable to the 
data presented here. Hatch et al. (2007) presents IFU obser- 
vations of six emission-line nebulae in cool-core clusters (se- 
lected from the optical ROSAT follow-up by Crawford et al. 
1999) with OASIS on the WHT (with a limited wavelength 
range centered around Ha). Edwards et al. (2009) presents 
IFU observations of nine CCGs in cooling and non-cooling 
clusters (also from Crawford et al. 1999) and within 50 kpc 
of the X-ray centre, with CMOS IFU and OASIS on WHT 
(most galaxies with wavelength ranges also limited to the 
Ha region, but three also included H/3). Emission line maps 
of morphology, kinematics, often line ratios, and the stel- 
lar continuum have been published for Abell 496 by Hatch 
et al. (2007), and for Abefl 2052 by Edwards et al. (2009). 
These studies concluded that no proposed heating mecha- 
nism reproduces all the emission-line properties within their 
observed wavelength within a single source. Thus, a single 
dominant mechanism may not apply to all CCG nebulae, 
and there may be a mixture of heating mechanisms acting 
within a single nebula (Wilman, Edge & Swinbank 2006; 
Hatch et al. 2007). Thus to get more information which 
will enable the dominant mechanism(s) to be identified, we 
observed more line ratios over a longer wavelength range 
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(around Ha and H/3 which also allowed extinction to be ap- 
propriately measured). We will be adding the detailed stel- 
lar population analysis, and place the derived information 
from the optical spectra in context with multiwavelength 
data over the full spectrum in a future paper, to explain the 
diverse nature of these galaxies. All the data for Abell 780 
and 1644 is new, and all four of the extinction maps are new. 
Also, the construction of ionisation diagrams and analysis of 
ionisation processes with shock models and pAGB models 
is also new. Our exposure times (improved from previous 
studies) are shown in Table [l] This data is also complimen- 
tary to the long-slit spectroscopy (on KECK and Magellan) 
along the Ha filaments of the four objects (studied here) by 
McDonald et al. (2012). 

We introduce the sample and detail of the data reduc- 
tions in Sections |2] and [3] We then derive the optical extinc- 
tion as well as the line strengths in Section |4]We proceed to 
discuss the four individual cases in Sectiori[5] and identify 
the mechanism producing the emission and the ionised gas 
in these four CCGs in Section [6] We summarise the findings 
of this paper in ^ We have used the following set of cos- 
mological parameters: Qrn — 0.3, Qa = 0.7, Ho = 70km s~^ 
Mpc-\ 



2 SAMPLE 

We have chosen our sample of active central cluster galaxies 
from the Ha imaging presented in McDonald et al. (2010), 
who in turn, selected their sample from White, Jones & 
Forman (1997). McDonald et al. (2010) enforced the cuts: 
6 < +35deg and 0.025 < z < 0.092, after which they se- 
lected 23 clusters to cover the full range of properties, from 
very rich clusters with high cooling rates to low-density clus- 
ters with small cooling flows. Their classical cooling rates 
range from 6.3 - 431 M© yr"^ which means that while cov- 
ering a large range in properties, their sample consisted of 
only cooling flow clusters. From their 23 cooling flow clus- 
ters, we selected all the clusters with clearly detected Ha in 
their centres (albeit filamentary, extended or nuclear emis- 
sion). In addition, all of these central galaxies have optical 
imaging, near-IR (2MASS) and UV (Galex data) available. 
Thereafter, we selected all the central galaxies with detailed 
X-ray (Chandra) data, as well as VLA 1.4 GHz fluxes, avail- 
able. This resulted in a sub-sample of 10 galaxies. We ob- 
served four of these galaxies with the CMOS IFU (as shown 
in Figure |2|. We merely chose the objects with the most 
auxiliary information available. This additional information 
will be added in the future paper (where the underlying stel- 
lar populations will be analysed in detail) and might help to 
constrain the ionisation mechanisms. 

The rest-wavelength range of the emission lines of in- 
terest is 4860-6731 A (H/3 to [Sn]A6731). The ratio of the 
forbidden [NH]A6584 to Ha line will depend on the metal- 
licity of the gas, the form of the ionising radiation, and the 
star formation rate. The relative strength of the [OHI]A5007 
and H/3 lines reveals further excitation mechanism and gas 
metallicity information. The role of AGN photoionisation is 
confined to the central 2 - 3" of active, massive nearby ellip- 
tical galaxies (Sarzi et al. 2006). Thus, IFU observations are 
ideal and will also allow us to study the 2D-distribution of 
the ionising radiation. In addition to the information from 



the emission lines, we are able to extract the underlying stel- 
lar absorption spectra using the improved GANDALF code 
(Sarzi et al. 2006). Thus, the kinematics and morphology of 
the hot ionised gas and stellar components can be correlated. 



3 OBSERVATIONS AND DATA REDUCTION 

The data were obtained with the CMOS IFU on the Gemini 
South telescope in semester 2011 A (February to July 2011). 
The GMOS-IFU in 1-slit mode was used and allowed us to 
map at least a 3kpc wide region in the centre of the target 
galaxies with a simultaneus coverage of the 4600-6800 A 
range in the target rest frame (using the B600 grating) with 
a single pointing. This resulted in a spectral resolution of 1.5 
A. This spectral resolution (81 km s~^) is poorer than that of 
Edwards et al. (2007, who had a much shorter wavelength 
range), and much higher than that of Hatch et al. (2006) 
(223 - 273 km s"^). 

The IFU field-of-view is 5 x 3.5 arcsec, and this area is 
divided into 500 lenslets (and another 250 lenslets offset for 
sky measurements). We obtained six exposures per galaxy, 
with the exception of MCG-02- 12-039 where we obtained 7 
exposures, resulting in a total of 12500 galaxy spectra. The 
targets and exposure times are shown in Table ^ Our in- 
tegration time is three times that of Edwards et al. (2009), 
and more than five times that of Hatch et al. (2007) with a 
much bigger instrument. In addition to the targets, the nec- 
essary bias, flat-fields, twilight flat-fields, and arcs frames at 
two different central wavelengths were also observed, as well 
as a spectrophotometric standard star for flux calibration. 
Two central wavelength settings were used to avoid losing 
crutial spectral information in the two CCD gaps. 

For more detail on the CMOS IFU data reduction pro- 
cess see Gerssen et al. (2006). The basic data reduction was 
done using the CMOS package in IRAF. The IFU sky-to- 
detector mapping was stored in the data array prior to data 
reduction. Several bias frames were averaged and subtracted 
directly from raw data, to correct the zero point for each 
pixel. Additional care was taken to avoid including raw bias 
frames in the mean bias frame that drifted measurably with 
time. Frames were mosaiced, and the overscan regions were 
trimmed. Flat-field and twilight fiat-field frames were used 
to correct for differences in sensitivity both between detector 
pixels and across the IFU field. The majority of the cosmic 
rays were rejected in the individual frames before sky sub- 
traction using the Gemini cosmic ray rejection routine. The 
remainder of the cosmic rays were eliminated using the LA- 
Cosmic routine (van Dokkum 2001) with an IRAF script 
that retained the multi-extension fits format for further re- 
ductions. The sets of 2D spectra were calibrated in wave- 
length using the arc lamp spectra for the two different cen- 
tral wavelength settings. The IFU elements were, thereafter, 
extracted from the raw data format to a format more con- 
venient for further processing. Sky emission lines and con- 
tinuum were removed by averaging the sky spectrum over a 
number of spatial pixels (from the offset sky fibers on the 
edge of the science field) to reduce the noise level, before 
subtracting it from all the spatial pixels. Thus the process 
adds little extra noise to the result since the observations 
were obtained in dark time (the variability of the sky region 
was minimal), and a number of 250 spaxels were averaged 
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(a) Abell 0496 - MCG-02-12-039 (b) Abell 0780 - PGC026269 (c) Abell 1644 - PGC044257 



(d) Abell 2052 - UGC09799 




(e) Abell 0496 - MCG-02-12-039 (f) Abell 0780 - PGC026269 



Abell 1644 - PGC044257 



(h) Abell 2052 - UGC09799 







(i) Abell 0496 - MCG-02-12-039 (j) Abell 0780 - PGC026269 (k) Abell 1644 - PGC044257 (1) Abell 2052 - UGC09799 

Figure 1. DSS images of the four targets (north above, and east to the left for all images). The upper plots show 8x8 arcmin field of 
views, and the middle plots show the targets with the 5 x 3.5 arcsec IFU field of view overlayed. The top of the IFU FOV is indicated 
with a blue arrow. The lower plots show continuum images made from the IFU cubes (width of 50 A at 6350 A), smoothed spatially 
with a Gaussian with width 3 spaxels (which corresponds to 0.3 arcsec) and using the Sauron colourmap. 



Object 


Cluster 


Redshift 


Linear scale 


^off 


Tx 


Classical cooling rates 


Spectrally determined 


Exposure Time 






z 


(kpc/arcsec) 


(Mpc) 


(keV) 


(Mo/yr-i) 


(Mo/yr-i) 


(seconds) 


MCG-02-12-039 


Abeh 0496 


0.0329 


0.654 


0.031 


4.8 


134 


1.5 


7 X 1800 


PGC026269 


Abeh 0780 


0.0539 


1.059 


0.015 


4.7 


222 


7.5 


6 X 1800 


PGC044257 


Abeh 1644 


0.0474 


0.935 


0.009 


5.1 


12 


3.2 


6 X 1800 


UGC09799 


Abeh 2052 


0.0345 


0.685 


0.038 


3.4 


94 


2.6 


6 X 1800 



Table 1. Galaxies observed with the Gemini South telescope. All four galaxies show extended Ho; emission (McDonald et al. 2010). The 
cluster X-ray temperature (Tx) and classical cooling rates rates (M) are from White et al. (1997). The spectrally determined cooling 
rates are from McDonald et al. (2010). The values for Ro// are from Edwards et al. (2007), with the exception of PGC044257 which is 
from Peres et al. (1998). 
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Object 


Cluster 


Rest wavelength range 

A 


Foreground extinction (mag) 

E{B-Y) galactic 


Average extinction (mag) 
measured E(B-Y)total 


Radio flux 
mJy 


MCG-02-12-039 


Abell 0496 


4648 - 7540 


0.140 


0.425 


121 


PGC026269 


Abell 0780 


4743 - 7693 


0.042 


0.210 


40800 


PGC044257 


Abell 1644 


4713 - 7646 


0.071 


0.195 


98 


UGC09799 


Abell 2051 


4655 - 7552 


0.037 


0.460 


5500 



Table 2. Further properties of the CCGs observed on Gemini South. Radio fluxes are from the NVSS. 



in the sky-subtraction process. Thus, the error contribution 
of the sky-subtraction process is ^fg^x the error on one 
sky spaxel. A spectrophotometric standard star (LTT4816) 
was used to correct the measured counts for the combined 
transmission of the instrument, telescope and atmosphere 
as a function of wavelength. We reduced the standard star 
observation with the same instrument configuration as the 
corresponding scientific data. A ID spectrum was extracted 
by adding the central spatial pixels from the standard star 
observation, and it was used to convert the measured counts 
from the galaxy spectra into fluxes with erg cm~^ s~^ 
units. 

The reduced 2D arrays were transformed back to a 
physical coordinate grid (x,y,X datacube) before scientific 
analysis, while also correcting for atmospheric dispersion. 
The latter (also called differential refraction) causes the posi- 
tion of a target within the IFU field to vary with wavelength. 
This correction was necessary as data were taken at differ- 
ent airmasses throughout the observing nights. The spatial 
offset as a function of wavelength was determined using an 
atmospheric refraction model (given the airmass, position 
angle and other parameters) from the SLALIB positional 
astronomy library. Each hexagonal spaxel (each spatial ele- 
ment) was 0.2 arcsec, and this was subsampled onto a rect- 
angular grid of 0.1 arcsec per spaxel when the separate expo- 
sures were combined. The exposures were combined (median 
averaged) using a centroid algorithm to calculate the shift- 
ing in X and and also shifting in A for the exposures at 
two different wavelength settings. No additional cosmic rays 
were visible after the LACosmic task was run on the indi- 
vidual cubes, so cosmic rays were not removed during this 
step. The cubes were also converted into RSS (row-stacked 
spectra) for further data reductions in IDL. Each spaxel was 
averaged with its eight neighbouring spaxels (improving the 
S/N by a factor of three), which is effectively smoothing over 
0.3 arcsecs - this is still slightly undersampled compared to 
the average seeing (^1 arcsecs), but only larger regions will 
be analysed further. 



(a) 6548.2 A 




(d) 6556.2 A 




(g) 6564.1 A 




(b) 6550.9 A 




(e) 6558.8 A 




(h) 6566.7 A 




(c) 6553.5 A 




(f) 6561.4 . 




(i) 6569.3 . 



Figure 2. Slices of width 1 A (smoothed spatially with Gaus- 
sian FWHM 1.5 spaxels and using a SAURON colourmap) of 
PGC026269 showing the wavelength and spatial scale over which 
the transition from [NII]A6548 to Ho: emission occurs, and how 
the morphology of Ho; changes. Although the continuum emission 
is usually smooth, the morphologies of the line emission are not 
uniform. 



4 LINE MEASUREMENTS AND INTERNAL 
EXTINCTION 

To accurately measure the emission-line fluxes of the CCG 
spectra, we use a combination of the ppxf (Cappellari & 
Emsellem 2004) and GANDALF (Gas AND Absorption Line 
Fitting algorithm; Sarzi et al. 2006) routine^ Gandalf ver- 

^ We make use of the corresponding ppxf and GANDALF IDL 
(Interactive Data Language) codes which can be retrieved at 
http: / www.leidenuniv.nl/sauron/. 



sion 1.5 was used as it enables a reddening correction to be 
performed, and it encorporates errors. This code treats the 
emission lines as additional Gaussian templates, and solves 
linearly at each step for their amplitudes and the optimal 
combination of stellar templates, which are convolved by 
the best stellar line-of-sight velocity distribution. The stel- 
lar continuum and emission lines are fitted simultaneously. 
All 985 stars of the MILES stellar library (Sanchez-Blazquez 
et al. 2006) were used as stellar templates to automatically 
include a-enhancement in the derived optimal template. The 
emission lines were masked while the optimal template was 
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derived. The Ha and [NII]A6583 lines were fitted first, and 
the kinematics of all the other lines were tied to these lines, 
following the procedure described in Sarzi et al. (2006). How- 
ever, in cases where the emission of the other lines were 
strong enough to measure velocity and velocity dispersion 
(as was mainly the case except for the extreme edges of the 
data cube), this was calculated independently as there is no 
a priori reason to expect the kinematics measured from all 
the lines to be the same (as they can originate in different 
regions). After the kinematics are fixed, a Gaussian tem- 
plate is constructed for each emission line at each iteration, 
and the best linear combination of both stellar and emission- 
line templates (with positive weights) is determined. This is 
done without assuming line ratios, except in the case of dou- 
blets where their relative strength is fixed by the ratio of the 
corresponding transition probabilities. We have adapted the 
gandalf code to apply it to the GMOS IFU cubes for a longer 
wavelength range. All 1617 spaxels were collapsed together 
to obtain a ID spectrum per cube, thereafter all 985 stars 
for the MILES library were used to create a global optimal 
template for each galaxy. This global optimal template (and 
the stars it consisted of - to account for varying weights over 
the spatial region) was then applied to all 1617 spectra per 
cube. The spectral types of the stars that the global optimal 
templates consisted of (from the MILES library) are shown 
in Table m 



Some ellipticals contain dust in the centre that can be 
patchy, uniform or filamentary (Laine et al. 2003). The long 
wavelength range of the spectra allows us to constrain the 
amount of reddening using the observed decrement of the 
Balmer lines, which can be set to have an intrinsic decre- 
ment consistent with the recombination theory by treating 
the lines as a multiplet. The physical constraints on the 
emission from the higher-order Balmer lines also ensures 
the strength of the corresponding absorption features is cor- 
rectly estimated. 



We used the dust models by Calzetti et al. (2000) to 
calculate the flux attenuation values at the desired wave- 
length for any given E{B — V) value (optional see below). 
The Balmer decrement assumes a case B recombination for 
a density of 100 cm~^ and a temperature of lO'^ K, result- 
ing in the predicted Ha/H/3 ratio of 2.86 (Osterbrock 1989). 
The code can adopt either a single dust component, affect- 
ing both the stellar continuum and the emission- line fluxes, 
or in addition include a second dust component that affects 
only the emission- line templates. We did not specify Galac- 
tic extinction (from the NED database, Schlegel, Finkbeiner 
& Davis 1998), but we measured a total diffuse component 
which gave the total extinction (hence including the fore- 
ground galactic extinction). This total extinction measured 
from the Balmer decrement are shown in Figure [S] and av- 
erages are noted in Table [5] (the galactic extinction noted in 
Table [2] was not subtracted). The extinction was smoothed 
over 0.3 arcseconds, and is only plotted where the velocity 
dispersion of the Ha line is less than 500 km s~^ (to avoid 
spaxels where the Ha line could not be separated from the 
[Nil] lines), and it is only plotted where the amplitude-to- 
noise (A/N) ratio of the Ha line is higher than 3 (as defined 



in Sarzi et al. 2006^ The A/N maps of the Ha measure- 
ments are shown in Section |5] which gives a direct indication 
of the errors on the extinction. 

The parameter E(B — V), i.e. the colour excess between 
4350 A and 5550 A, for the galactic extinction for each of the 
four galaxies was taken from the NED database (Schlegel et 
a. 1998), and ranged between 0.037 and 0.140 mag (Table|2|. 
The parameter Rv, i.e. the ratio of the absolute extinction 
at 5550 A (Av) to the colour excess E{B — V), was taken as 
3.1 for the interstellar medium (Cardelli, Clayton & Mathis 
1989). The total extinction is given by: 



E{B - V)total = 



2.177 
-0.37Rv 



X (lo^ 



Hao 



log 



Haobs 
Hpobs 



The theoretical Ha/H/3 flux ratio of 2.86 may not be the 
ideal value to use for Seyfert-type galaxies, but the actual 
value is debated. It is often assumed that the Ha emission in 
these systems is enhanced due to collisional processes, and 
several authors use a value Rv of 3.1 (Gaskell & Ferland 
1984; Osterbrock & Ferland 2006), although other values 
have also been determined (Binette et al. 1990 calculate a 
value of 3.4). The total extinction maps are presented in 
Figure [3] and shows mostly very low extinction, but some 
morphological features can be seen in MCG-02- 12-039 and 
UGC09799. Particularly high values of E{B — V) internal can 
be seen in UGC09799. The galactic extinction of PGC044257 
is E{B — V)gai = 0.071 mag, and from long-slit spectra, 
Crawford et al. (1999) derived the total extinction as 0.46 
to 0.63 mag. This agrees with the extinction we derived for 
the very centre of the galaxy in Figure [3] but on average 
our spatially resolved extinction is slightly lower. The galac- 
tic extinction E{B - V)gai of UGC09799 is 0.037 mag, and 
Crawford et al. (1999) derived an integrated internal ex- 
tinction of E(B — V)internai of 0.22 mag for the centre of 
this galaxy. This corresponds very well to what we derived 
and plotted in Figure |3] although we find that some regions 
show much higher internal extinction. The values of extinc- 
tion determined here may be slightly overestimated due to 
the choice of intrinsic Ha/H/3 flux used. 

Figure|2]show slices with width 1 A (smoothed spatially 
with Gaussian FWHM 1.5 spaxels and using the Sauron 
colourmap) of PGC026269 showing the wavelength and spa- 
tial scale over which the morphology of Ha changes. The 
morphology of the nuclear region changes rapidly across the 
emission lines. This is illustrated in the figure showing a se- 
ries of monochromatic slices. Although the continuum emis- 
sion is usually smooth, the morphologies of the line emission 
are not uniform. 

We are able to measure Ha line fluxes with 3 A/N 
accuracy at around 5 x 10~^^ erg cm~^ s~^. After cor- 
recting for extinction, we are able to measure the fol- 
lowing lines within our wavelength range: [ArIV]A4740 
A; H/3A4861 A; [OIII]AA4958,5007 A; [NI]AA5198,5200 A; 
HeIA5876 A; [OI]AA6300,6364 A; [NII]AA6548,6583 A; Ha 
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at 6563 A; [SII]AA6716,6731 A (see Figures [e] to |2q| and 
for the A/N ratios of the individual emission lines). Ha and 
[NII]AA6548,6583 A has already been measured for MCG- 
02-12-039 and UGC09799, and H/3 and [OIII] AA4958,5007 



The A/N is related to the S/N: EW ■ 



S 



A/ATx \/27rcr 
SjN 



where 



EW is the equivalent width of the line, and a is the line width. 
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Object Spectral types in the global optimal template 



MCG-02-12-039 AO la, GO, Gl lb, G5, G8 III, KO, KO III, K IIvw 

PGC026269 AO la, KO Ibpvar, M4 III 

PGC044257 B3 III, A5, F3 V, F6, F7 V, GO Vw, Gl lb, G5, G5 V, G8 III/IVw, Kl III 
UGC09799 B8 lb, A5, F3 III, GO Vw, Gl lb, KO III 



Table 3. Properties of the underlying stellar populations. 

A has also been measured for UGC09799. All the other line 
measurements are new. 
.pdf 



5 INDIVIDUAL GALAXIES 

The continuum emission of the galaxies are shown in Figure 
|2] Figure [4] [8] [12] and [l6] shows a random spectrum in the 
central region of MCG-02-12-039, PGC026269, PGC044257 
and UGC09799 respectively. The best-fitting stellar tem- 
plate, Gaussians at the emission lines, best-fitting stellar 
templates with the emission lines subtracted and the rel- 
ative flux for the measured emission lines are also shown. 
The absorption-extracted, dereddened maps of the Ha, 
[NII]A6583/Ha, and [OIII] A5007/H/3 emission are shown in 
Figure [BT] 

For comparison purposes, the kinematics were extracted 
from the IFU image along the same slit position as the long- 
slit data in Loubser et al. (2008). This comparison for MCG- 
02-12-039, PGC026269 and PGC044257 is showed in Section 
[A] and the data points compare satisfactory. 

The absorption-extracte d , de reddened maps of the Ha 
emission are shown in Figure |B1| in units of 10~^^ erg cm~^ 
s~^. We also show the useful ratios [NII]A6583/Ha, and 
[OIII]A5007/H/3. Figures [6] to [20] shows the A/N ratios of 
the Ha, [NII]A6583 A; [SII]AA6716,6731 A; [OIII]A5007 A; 
[OI]A6300 A lines. As mentioned in the Introduction, MCG- 
02-12-039 and UGC09799 already has previous IFU observa- 
tions which is compared (in terms of S/N and lines measur- 
able) in the beginning of Section 3 and the last paragraph 
of Section 4. The stellar kinematics in Figures [5] [9] [13] and 
|18| were measured using all the absorption lines within our 
wavelength range. 

For a sample of ^50 elliptical galaxies, Sarzi et al. 
(2006) find ionised gas velocities (estimated using the 
[OIII]A5007 line) between -250 and 250 km s"\ and gas 
velocity dispersion as high as 250 km s~^. Using the same 
line, we find gas velocities from ±100 to ±350 km s~^, and 
line widths from 200 to 420 km s~^. Using Ha, we find gas 
velocities from ±125 to ±350 km s~^, and line widths from 
200 to 400 km s"\ 



5.1 MCG-02-12-039 

The central galaxy (MCG-02-12-039) of Abell 0496 is a fairly 
weak line emitter (Fabian et al. 1981; Hu, Cowie & Wang 
1985) and is host to the compact radio source MSH 04-112 
(Markovic, Owen k Eilek 2004). Abell 0496 is a relaxed 
cluster with a cool core with a central metal abundance en- 
hancement (Tamura et al. 2001). At a redshift of 0.0329, the 
galaxy has a linear scale of 0.656 kpc arcsec"^ (from Hatch 



et al. 2007). This cluster has an interesting Ha morphology 
(McDonald et al. 2010). There are at least 5 distinct fila- 
ments, with various shapes and directions. The two longest 
filaments run parallel to each other for ^12 kpc (18.3 arc- 
sec, whereas our observations only cover 3.5 by 5 arcsecs in 
the centre). 

Figure |B1| shows that the peak Ha emission is also 
where the dust (extinction) is the highest in Figure [s] This 
is similar to the finding of Hatch et al. (2007), where they 
found that the peak Ha ± [Nil] emission corresponds to the 
area where the three dust lanes meet in the galaxy centre. 
We confirm the observation by Hatch et al. (2007), that the 
line emission follows the general path of the dust (extinc- 
tion) features, but is not as filamentary. 

IFU observations of this galaxy was also taken by 
Hatch et al. (2007), but they only measured the Ha and 
[NII]AA6548, 6583 A lines. Their spectral resolution is much 
poorer at 223 - 273 km s~^ than ours at 81 km s~^. Our 
S/N is also more than five times that of Hatch et al. (2007), 
observed with a much bigger instrument. We plot Ha, [Nil], 
[SII], [OIII] and [01] velocities, line width and A/N in Fig- 
ure [6] and [7] The stellar component displays a elongated 
morphology in the kinematics plotted in Figure [5] The gas 
components also appear elongated (kinematics plotted in 
Figure |6|, but to a lesser extend. It suggests that the stars 
and gas are kinematically decoupled. Comparison of the gas 
kinematics (Figures [6] and [7| with the stellar kinematics pre- 
sented in Figure[5] and Figure [AT] as well as by Fisher, Illing- 
worth & Franx (1995) confirms the Hatch et al. (2007) ob- 
servations that the two components are disconnected. Com- 
parison of the plots in Figure [6] showed that, to the degree 
that our spatial resolution reveals, it appears that all the op- 
tical forbidden and hydrogen recombination lines originate 
in the same gas. 

We find a maximum emission line width of 400 km s~^ 
as shown in Figures [6] and [7] This is again similar to the find- 
ing of Hatch et al. (2007). They found a maximum linewidth 
of 600 km s~^ in the dust-free central region to the north- 
east of the galaxy centre, and the rest of the nebula to have 
a linewidth of 100 - 250 km s~\ 

The stellar line-of-sight velocity reveals that some parts 
of the nucleas are blueshifted by ~125 km s~^ and some 
parts are redshifted by the same amount. Hatch et al. (2007) 
found the southern part of the galaxy to be blueshifted by 
-200 km s~^ whilst the northern section is marginally red- 
shifted up to ±150 km s~^. No clear kinematic pattern is 
associated with the dust structures. 

The peak- to-peak gas velocity of 250 km s~^ is fairly 
low compared to the other CCGs, and the stellar component 
of Abell 0496 has a mean rotation of 59 km s~^ (Paper 1). 



8 Louhser et al 




9 



0.012 



0.010 - 



c 0.008 - 

(D 



0.006 ■ 



CD 



0.004 ■ 



0.002 



0.000 I , .1 , , , I , l\ , I , , A 

1200 1400 1600 1800 2000 
Pixels 



0.012 



0.010 - 



c 0.008 

(D 



0.006 



0.004 



0.002 



0.000 



4600 4800 5000 5200 5400 
Pixels 



Figure 4. A random spectrum in the central region of MCG-02-12-039. The red hne indicates the best-fitting stellar template and 
Gaussians at the emission lines, and the green line indicates the best-fitting stellar templates with the emission lines subtracted. The 
blue line indicates the relative fiux for the measured emission lines. The lines are (from left to right) in the left plot: [ArlV], H/3, and the 
[OIII] doublet, and in the right plot: [OI] doublet, [Nil] doublet, and Ho;. The Ho; A/N of this pixel was 15. 
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Figure 8. A random spectrum in the central region of PGC026269. The red hne indicates the best-fitting stellar template and Gaussians 
at the emission lines, and the green line indicates the best-fitting stellar templates with the emission lines subtracted. The blue line 
indicates the relative fiux for the measured emission lines. The lines are (from left to right) in the left plot: H/3, the [OIII] doublet, and 
[NI], and in the right plot: [OI] doublet, [Nil], Ho;, [Nil], and the [SII] doublet. The sharp feature to the left of the first [Nil] line is one 
of the two CCD chip gaps and was masked during the pPXF and GANDALF fitting processes. 
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5.2 PGC026269 

This galaxy coincides with the well-known luminous radio 
counterpart Hydra A. Abell 0780 is a poor cluster with an 
associated cooling flow nebulae. IFU observations of this 
galaxy has not been presented in the literature before this 
study. 

The stellar line-of-sight velocity (from the IFU data) is 
blue- and redshifted by ^175 km s~^ which is slightly higher 
than that derived from long slit spectroscopy (51 ± 20 km 
s~^ Paper 1). 

We plot the kinematics of the stellar and gaseous com- 
ponents in Figures |9] and |10| Both the stellar and gaseous 
components show clear rotation in Figures [9] and (To] al- 
though it also seems to be decoupled. This system bears 
resemblance to the BCG NCG 3311 in Abell 1060, studied 
by Edwards et al. (2009) which also showed striking rota- 
tion. The rotation of the warm gas agrees with the long-slit 
observations by McDonald et al. (2012). 

We plot Ha, [NH], [SH], [OHI] and [01] velocities, line 
width and A/N in Figure [To] and 11 Comparison of the 
plots in Figure [6] showed that, to the degree that our spatial 
resolution reveals, it appears that all the optical forbidden 
and hydrogen recombination lines originate in the same gas. 

Wise et al. (2007) present a summary of the X-ray and 
radio properties of this cluster, showing the excellent cor- 
relation between the radio jets and the X-ray cavities. The 
arcing Ha filament that McDonald et al. (2010) detect north 
of the CCG (on a larger scale as this study) appears to be 
spatially correlated with the radio jet. Our Figure 24 (d - 



f) show that our optical observations confirms that the pho- 
toionisation is caused by an AGN (the green and red curves) , 
and are therefore consistent with the X-ray and radio studies 
of this system on larger scales. 



5.3 PGC044257 

This CCG is one of the rare cases (occurs in ^ 3 per cent of 
CCGs; Hamer et al. 2012) where the CCG is slightly offset 
from the optical line emission (Johnson et al. 2010). The 
rarity of such offsets points to a large event in cluster evo- 
lution, a major cluster merger or possibly a powerful AGN 
outburst. Whatever the reason for the separation the gas 
cooling at the X-ray peak will continue and cooled gas will 
be deposited away from the CCG (Hamer et al. 2012). 

IFU observations of this galaxy has not been presented 
in the literature before this study. The stellar line-of-sight 
velocity is blue- and redshifted by ^225 km s~^ which is 
higher than derived from long slit spectroscopy (20 ± 16 km 
s~^ Paper 1). 

We plot the kinematics of the stellar and gaseous com- 
ponents in Figures [13] and [M] Both the stellar and gaseous 



components show morphologies (Figures 13 and 14) that 



are elongated and aligned. However the recession velocity of 
the gas is ^ 200 km s~^, which is higher than that of the 
stars, suggesting that the gaseous and stellar components 
are decoupled. The rotation of the warm gas agrees with 
the long-slit observations by McDonald et al. (2012). 

We plot Ha, [Nil], [SH], [OHI] and [01] velocities, line 
width and A/N in Figure Il4] and 15 To the degree that 



forbidden and hydrogen recombination lines originate in the 
same gas. 

We find the Ha flux to be uniform in Figure |B1| sug- 
gesting that all of the gas had the same origin, but the Ha 
flux is also quite low. If Figure [BT] is compared to Figure [M] 
then it can be seen that some structure is visible where the 
A/N of Ha is the highest. The Ha appear to be quite core- 
dominated. McDonald et al. (2010) found that the direction 
of the Ha filaments (from narrow band imaging) correlate 
with the position of nearby galaxies. The filaments in Abell 
1644 also appear to be star- forming (McDonald et al. 2011). 
The star forming filaments that McDonald et al. (2011) de- 
tected are on a bigger scale than the region we probe in this 
study. Our results show that the very centre of the system 
display properties of a LINER (see Section [6]), and therefore 
both mechanisms may be part of this system at different 
scales. 

5.4 UGC09799 

IFU observations of this galaxy (around Ha and H/3) was 
previously presented by Edwards et al. (2009). The obser- 
vations presented in this study are improved in that it was 
observed with triple the integration time of the previous ob- 
servations. The morphology of the continuum emission of 
this galaxy (see Figure [2]) is centrally concentrated and con- 
densed, as was also found by Edwards et al. (2009). 

Just fitting one single Gaussian per emission line re- 



sulted in very poor fits (as shown in Figure 17). We also 
tried fitting Voigt profiles to the [OIII] lines, for example, as 
well as slightly offset blue and red velocity wings. The best 
fits were consistently achieved using 2 Gaussians for each 
individual hne in the [OIII], [01], [Nil] and [SH] doublets, 
of which one of the Gaussians were broader than the other. 
The H/3 lines (and therefore also Ha) however required no 
additional Gaussian to achieve a good fit. Edwards et al. 
(2009) could not detect H/3 emission above the la level for 
this galaxy. We do detect H/3 above this level and plot the 



flux ratios in Figure Bl Figure 16 in Edwards et al. (2009) 
can be compared with Figures ^] and |B1| here - both show 
smooth centrally condensed emission. 

We plot the kinematics of the stellar and gaseous com- 



ponents in Figures [T8| and [T9| We plot Ha, [Nil], [SH], [OIII] 
and [01] velocities, line width and A/N in Figure and [20) 
To the degree that our spatial resolution reveals, it appears 
that all the optical forbidden and hydrogen recombination 
lines originate in the same gas. Our kinematic analysis can 
be compared to that of Edwards et al. (2009), who detected 
a gradient in the velocity but could not differentiate be- 
tween a rotation or outflow. Our stellar kinematics is shown 
in Figure [18] however we do not detect rotational kinemat- 
ics. Edwards et al. (2009) also found Ha to have a velocity 
range of -250 km s~^ to +150 km s~^. Our Ha kinemat- 
ics are shown in and show also show a range (^ 400 km 



our spatial resolution reveals, it appears that all the optical 



s~^) of velocities for Ha. The gaseous components show ro- 
tation in Figure [19] even though no rotation is apparent in 
the stellar components (Figure [iS]). 

This galaxy has patchy dust in the centre (Laine et 
al. 2003, also shown in Figure |3|, and Hicks & Mushotzky 
(2005) have deduced star formation from the excess UV- 
IR emission. The FUV/Ha ratio suggests heating by fast 
shocks or some other source of hard ionisation (e.g., cosmic 
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Figure 12. A random spectrum in the central region of PGC044257. The red hne indicates the best-fitting stellar template and Gaussians 
at the emission lines, and the green line indicates the best-fitting stellar templates with the emission lines subtracted. The blue line 
indicates the relative fiux for the measured emission lines. The lines are (from left to right) in the left plot: H/5, and the [OIII] doublet, 
and in the right plot: [OI] doublet, [Nil], Ha, [Nil], and the [SII] doublet. 
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Figure 16. A random spectrum in the central region of UGC09799. The red hne indicates the best-fitting stehar template and Gaussians 
at the emission lines, and the green line indicates the best-fitting stellar templates with the emission lines subtracted. The blue line 
indicates the relative fiux for the measured emission lines. The lines are (from left to right) in the left plot: H/5, and the [OIII] doublet, 
and in the right plot: [OI] doublet, [Nil], Ho;, [Nil], and the [SII] doublet. The sharp feature between the [Nil] and [SII] lines is one of 
the two CCD chip gaps and was masked during the pPXF and GANDALF fitting processes. 
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(a) Single Gaussians 



(b) Two Gaussians 



Figure 17. UGC09799 showing single Gaussian fits and double Gaussian fits (one narrow and one broad) to the [OI] doublet, [Nil], Ho;, 
[Nil], and the [SII] doublet. 



rays, AGN). The filaments of Abell 2052 have ratios which 
are consistent with this picture (McDonald et al. 2011). The 
highest extinction (in Figure [3| does not coincide with the 
Ha peak (Appendix Bl), and is slightly ofl[ centre. This does 
not necessarily contrast with the results quoted by Laine 
et al. (2003) that the dust is patchy in the centre, as the 
centre of the observations might not necessarily be exactly 
the same as our IFU placement (see Figure 1). Figure 1 also 
shows that our continuum emission is smooth, similar to the 
conclusion reached by Edwards et al. (2009). 



Venturi, Dallacasa & Stefanachi (2004) found a parsec- 
scale bipolar radio source, and Chandra X-ray emission show 
two bubbles in the ICM on a larger scale (Blanton, Sarazin & 
McNamara 2003). Our results are consistent with the central 
emission being that of a LINER (see Section [gJ. 



6 DISCUSSION - lONISATION MECHANISMS 

Various mechanisms have been proposed as sources of exci- 
tation in the nebulae, including photoionisation by radiation 
from the AGN, cluster X-rays, or hot stars, collisional heat- 
ing by high-energy particles, shocks or cloud-cloud collisions, 
and conduction of heat from the X-ray corona (Peterson et 
al. 2003). However, none so far satisfactorily describes all 
the characteristic spectra, energetics, and kinematics of the 
extended emission-line regions (Wilman et al. 2006; Hatch 
et al. 2007). 

We reiterate that the FOV of the observations are only 
3.5 X 5 arcsec (corresponding to the central few parsec - see 
Table [T]). These results only look at the very heart of these 
sources. Therefore, it is really only the emission mechanisms 
for the gas at the very center of these systems that is being 
studied here. 

The forbidden lines such as [NII]A6583 result from the 
excitation of N^ through collisions with electrons liberated 
through photoionisation. The Ha emission results from the 
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recombination of the hydrogen ion. [NII]A6583 flux depends 
on the abundance, the strength of the radiation field, 
and the form of the radiation field: a harder ionising source 
will produce a greater fiux. The Ha fiux also depends on 
the strength of the radiation field. Therefore, the ratio will 
depend on the metallicity of the gas and the form of the ion- 
ising radiation. The form of the ionising radiation and/or the 
gas metallicity are not uniform but must vary within each 
galaxy and between the whole sample. Different excitation 
mechanisms may act in different regions (Hatch et al. 2007, 
Edwards et al. 2009), or the Ha emission might be disturbed 
by the presence of companion galaxies (Wilman et al. 2006). 

A commonly used method to distinguish between the 
sources of ionisation uses the emission-line diagrams pi- 
oneered by Baldwin, Phillips & Terlevich (1981, here- 
after BPT diagram) which separate the two major ori- 
gins of emission: star formation and AGN. The dia- 
grams use pairs of emission line ratios, of which the most 
commonly used is [OHI] A5007/H/3 against [NH]A6584/Ha. 
Extinction-corrected emission-line measurements were used 
(even though the BPT diagram is almost insensitive to red- 



dening). We plot the BPT diagrams in Figures 22 and 23 
using the following criteria: 

Kewley criteria: Galaxies above this line are AGN. Kew- 
ley et al. (2001) used a combination of photoionisation and 
stellar population synthesis models to place a theoretical up- 
per limit on the location of star forming galaxies on the BPT 
diagram, log(^^^ = .61 /\og{[N 1 1]/ Ha) - 0.47+1.19). 

Kauffmann et al. (2003) criteria: Galaxies below this 
line are purely star forming. Between Kewley and Kauff- 
mann criteria is composite objects. Kauffmann et al. (2003) 
revised the criteria as follows: A galaxy is defined to be an 
AGN ff log(^^^ > OM/\og{[NII]/Ha) - 0.05 + 1.3). 

LINERS-Seyffert line (Schawinsky et al. 2007): 
This line distinguishes LINERS from Seyffert galaxies 
log(I^^) = 1.051og(I|^^l) + 0.45. 

SDSS emission-line galaxies occupy a well-defined re- 
gion shaped like the wings of a seagull (Stasinska et al. 2008). 
Although the exact location of the line dividing the star 
forming and active galactic nuclei galaxies is still controver- 
sial (Kewley et al 2001; Kauffmann et al 2003; Stasinska et 
al 2006). 

All galaxies show important LINER emission, but that 
at least one has significant Seyfert emission areas, and at 
least one other has significant HII like emission line ratios 
for many pixels as shown in Figures [22] and |23| This is in 
agreement with the long-slit data of these sources plotted 
on BPT diagrams in figure 3 of McDonald et al. (2012). 

However, there is a debate about the ionisation mecha- 
nism in LINERS (Low-Ionisation Nuclear Emission-line Re- 
gion). The most viable excitation mechanisms are: a low 
accretion-rate AGN (Kewley et al. 2006), photoionisation by 
old post-asymptotic giant branch (pAGB) stars (Stasinska 
et al. 2008), and fast shocks (Dopita & Sutherland 1995). 
Sarzi et al. (2010) investigate the ionising sources for the 
gas in elliptical galaxies based on SAURON intergral-field 
spectroscopy whose spectra are limited to a relatively nar- 
row wavelength range. They conclude that pAGB stars are 
the main source of ionisation. In contrast, Annibali et al. 
(2010) analyse long-slit spectra of 65 ellipticals and claim 
that their nuclear line emission can be explained by excita- 
tion from the hard ionising continuum from an AGN and/or 



fast shocks. However, they can not completely rule out a con- 
tribution from pAGB stars at large radii. Voit & Donahue 
(1997) suggested that sources of supplementary heating pro- 
duce the LINER-like properties of the spectra, though not 
necessarily through the same mechanism in all systems. 

Annibali et al. (2010) found that from the centre out- 
ward, galaxies move left and down in the BPT diagram for 
their study of 65 early-type galaxies. Thus, the hardness 
of the ionising continuum decreases with galactocentric dis- 
tance (up to half the half-light radius in the Annibali et al. 
(2010) sample). Figures 22 and 23 also show the fiux ratios 



as function of distance from the galaxy centre. The red cir- 
cles indicate the central 0.5 x 0.5 arcsec of the galaxy, the 
yellow circles 1.0 x 1.0 arcsec, and the blue circles the full 
3.5 X 5.0 arcsec. The centre of the galaxy was determined as 
the luminosity peak in the continuum images in Figure [2] In 
our case, the hardness of the ionising continuum stay mostly 
uniform with galactocentric distance over out limited spatial 
extend. PGC044257 show an interesting core separation of 
the emission in the very centre of the galaxy in Figure [23^. 

There is also the debate about which galaxies are LIN- 
ERS (AGN) and which have just LINER-like emission (non- 
AGN) . The division is usually made by looking at the extend 
of the LINER signature: core dominated mean true LINER 
and diffuse means LINER-like. The arising problem is that 
the scale of the SDSS fibers is already too large to make 
the distinction for most galaxies. For the nearby galaxies, 
the centroiding of the fibres is not accurate enough to be 
sure if the measured spectra cover the core of the galaxy. 
We have core dominated LINER emission for at least three 
out of the four galaxies. Additional confirmation comes from 
the fact that at least one of our galaxies (PGC026269) is 
a strong radio source (AGN, see [2]). Similarly, our find- 
ings for UGC09799 agree with Edwards et al. (2009) who 
found line ratios consistent with Seyffert or LINER activ- 
ity in most of their central spaxels for UGC09799. Edwards 
(2009) also found that the Seyffert signature dominated the 
central spaxels of the CCG (UGC09799). The Ha emission 
surrounding the BCG in this cluster is coincident with radio- 
blown bubbles in the central region of the cluster. These bub- 
bles to the north and south of the cluster core are filled with 
radio emission, which likely originated from the AGN within 
the CCG (Blanton et al. 2003). Since the Ha emission seen 
by McDonald et al. (2010) is primarily along the edges of the 
northern bubble, they suspect that shocks may be respon- 
sible for the heating in this case. However, for PGC044257, 
McDonald et al. (2010) find very little evidence for an AGN 
in the cluster hosting the galaxy, in terms of the radio power. 
X-ray morphology and hard X-ray fiux. Two of the galaxies 
in the current study have detected X-ray point sources in 
the CCG (PGC026269 and UGC09799) with large cavities 
in their X-ray haloes, suggesting that the AGN is infiuencing 
the surrounding medium. 

We investigate this further by plotting the other 
BPT diagrams, on a spaxel-by-spaxel basis, for all 
four galaxies (shown in Figure |24|: [OIII] A5007/H/3 vs 
[NII]A6584/Ha, [OIII]A5007/H/3 vs [SII]AA6717,6731/Ha, 
and [OIII]A5007/H/3 vs [OI]A6300/Ha. One BCG 
PGC026269 show several HII pixels in all three BPT 
diagrams. Figure [22] show that these HII pixels occur 
throughout the centre of the galaxy. 

On these plots, we compare our observations to the pho- 
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toionisation models for pAGB stars with 7a—7aq (Binette et 
al. 1994). These models are consistent with most of our ob- 
servations. A model with Z=l/3 Z© will be shifted towards 
lower values on the x-axis of the three BPTs. The pAGB sce- 
nario has recently been revisited by Stasmska et al. (2008), 
whose extensive grid of photoionisation models (see their 
Figure 5) cover most of the regions occupied by our spa- 
tially resolved measurements. 

Three other grids of ionisation models are overplotted 
on the BPT diagrams (Figure [241. The plotted AGN pho- 



toionisation models (Groves, Dopita & Sutherland 2004) 
have an electron density, Ue = 100 cm~^, metallicities of so- 
lar, 7a—7aq (red grids), and twice solar (green grids), a range 
of ionisation parameter (-3.6 < log [/ < 0.0) and a power- 
law ionising spectrum with spectral index a = -2, -1.4, and 
-1.2. A harder ionising continuum, with a — -1.2, boosts 
[SII]A, A6717,6731 and [OI]A6300 relative to Ha. 

We also compared our results with shock models (Allen 
et al. 2008, purple grids). In Figure 24 we plot the grids 
with Z=Z0, preshock densities of 100 cm""^, shock veloci- 
ties of 100, 500 and 1000 km s~^, and preshock magnetic 
fields of B=l, 5 and 10 /xG. Shock models with a range of 
magnetic field strengths (B=l,5 and 10 /iG) match our ob- 
servations. Interstellar magnetic fields of B ^ 1 - 10 /iG are 
typical of what is observed in elliptical galaxies (Mathews & 
Brighenti 1997). Overall, shock models reproduce the ma- 
jority of our data in the three emission- line ratio diagrams. 
The shock grids with lower metallicity (e.g. LMC and SMC 
metallicities) are not consistent with our measurements |^ 

As shown in Table [2] - we have very weak as well as 
strong radio fiuxes in our small sample. We therefore believe 
that we are not particularly prone to biases such as the fact 
that a priori choice of galaxies with strong radio fiuxes will 
result in finding a sample where high [Nil] /Ha would be 
more common than in an optically, or Ha selected sample 
only - where starforming high Ha/ [Nil] sources might be 
more abundant. 



7 CONCLUSION 

We present detailed integral field unit (IFU) observations of 
the central few kiloparsecs of the ionised nebulae surround- 
ing four active CCGs in cooling fiow clusters (Abell 0496, 
0780, 1644 and 2052). Our sample consists of CCGs with 
Ha filaments. We observed the detailed optical emission-line 
(and simultaneous absorption line) data over a broad wave- 
length range to probe the dominant ionisation processes, 
excitation sources, morphology and kinematics of the hot 
gas (as well as the morphology and kinematics of the stars). 
Two of the four sources have not been observed with IFU 
data before (Abell 0780 and 1644), and for the other two 
sources we observed with significantly improved integration 
times (and number of lines) than previous studies (Hatch et 
al. 2007, Edwards et al. 2009). This will help form a com- 
plete view of the different phases (hot gas and stars) and how 
they interact in the processes of star formation and feedback 
detected in central galaxies in cooling fiow clusters, as well 
as the infiuence of the host cluster. 



^ We downloaded the shock and AGN grids from the web page 
ht t p : / / www . st rw . leidenuni v . nl / ^ brent /it er a . ht ml [ 



The total extinction maps are presented in Figure |3] 
and shows extinction which agrees well with previously de- 
rived long-slit values (where available). From long-slit spec- 
tra, Crawford et al. (1999) derived the total extinction of 
PGC044257 as 0.46 to 0.63 mag. This agrees with the ex- 
tinction we derived for the very centre of the galaxy in Figure 
|3] but on average our spatially resolved extinction is slightly 
lower (0.195 mag). Crawford et al. (1999) derived an inte- 
grated internal extinction of E{B — V)internai of 0.22 mag 
for the centre of UGC09799. This corresponds very well to 
what we derived and plotted in Figure [3] although we find 
that some regions show much higher internal extinction (on 
average 0.42 mag). 

We derive a range of different kinematic properties, 
given the small sample size. For Abell 0496 and 0780, we 
find that the stars and gas are kinematically decoupled, and 
in the case of Abell 1644 we find that these components are 
aligned. For Abell 2052, we find that the gaseous compo- 
nents show rotation even though no rotation is apparent in 
the stellar components. To the degree that our spatial reso- 
lution reveals, it appears that all the optical forbidden and 
hydrogen recombination lines originate in the same gas for 
all the galaxies. 

All galaxies show important LINER emission, but that 
at least one has significant Seyfert emission areas, and at 
least one other has significant HII like emission line ratios 
for many pixels (consistent with the long-slit observations of 
McDonald et al. (2012)). We also show that the hardness of 
the ionising continuum does not decrease with galactocen- 
tric distance (except for PGC044257 which show an inter- 
esting core separation of the emission in the very centre of 
the galaxy in Figure [23^). The radial profiles of diagnostic 
line ratios, [OIII]A5007/H/3 and [NII]A6584/Ha, show that 
they are roughly constant with radius for three of the four 
galaxies (all except PGC044257). This indicates that the 
dominant ionising source is not confined to the nuclear re- 
gion in the two objects and that the ionised gas properties 
are homogeneous in the emission-line regions across each 
galaxy. 

Overall, it remains difficult to dissentangle the domi- 
nant photoionisation mechanisms, even with more line mea- 
surements. The AGN photoionisation models (with higher 
metallicity) are the best able to reproduce our spatially- 
resolved line ratios the in all of the three BPT diagrams 
simultaneously of most objects, even though shock models 
and pAGB stars can not be conclusively eliminated. We also 
do not see extended [OI]A6300 emission, following the mor- 
phology of the strong emission lines (e.g. Parage et al. 2010), 
therefore, it is unlikely that shock excitation is the domi- 
nant ionising source in the galaxies of this limited sample. 
In addition, mult i wavelength observations (as discussed in 
the previous section) favour the AGN photoionisation mech- 
anism, especially in the case of PGC026269 and UGC09799. 
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Figure 22. BPT diagram for MCG-02-12-039 and PGC026269. The red circles indicate the central 0.5 x 0.5 arcsec of the galaxy, the 
yellow circles 1.0 x 1.0 arcsec, and the blue circles the full 3.5 x 5.0 arcsec. The centre of the galaxy was determined as the luminosity 
peak in the continuum images in Figure [2] The black solid curve is the theoretical maximum starburst model from Kewley et al. (2001), 
devised to isolate objects whose emission line ratios can be accounted for by the photoionisation by massive stars (below and to the 
left of the curve) from those where some other source of ionisation is required. The black-dotted curve in the diagram represent the 
Seyfert-LINER dividing line from Kewley et al. (2006) and transposed to the [Nil] A6584/Hq; diagram by Schawinski et al. (2007). 
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slit data in Paper 1. This comparison is showed in Figures 
|Al| to [A3l and the data points compare satisfactory. We em- 
phasize that the long-slit data were binned to a constant 
(higher) S/N, where the IFU data are three binned spaxels 
(0.3") along the same alignment. 



APPENDIX B: LINE MAPS 

This paper has been typeset from a TgX/ file prepared 

by the author. 



APPENDIX A: STELLAR KINEMATICS 
COMPARISON 

For comparison purposes, the kinematics were extracted 
from the IFU image along the same slit position as the long- 
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Figure Al. MCG-02-12-039: Kinematic data (velocity and velocity dispersion) extracted from the IFU data (empty symbols), compared 
with long-slit data along the same axis from Paper 1. 
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Figure A2. PGC026269: Kinematic data (velocity and velocity dispersion) extracted from the IFU data (empty symbols), compared 
with long-slit data along the same axis from Paper 1. 
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Figure A3. PGC044257: Kinematic data (velocity and velocity dispersion) extracted from the IFU data (empty symbols), compared 
with long-slit data along the same axis from Paper 1. 
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(d) PGC026269: Ha 



(e) PGC026269: [Nil] A6583/Ha 



(f) PGC026269: [OIII] A5007/H^ 



